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On the Molecular-Aerosol Scattering Coupling in 
Remote Sensing of Aerosol From Space 


V. V. Rozanov and A. A. Kokhanovsky 


Abstract—This work is devoted to studies of the influence of the 
molecular-aerosol coupling on the scalar approximation-based 
aerosol satellite remote sensing procedures in the ultravio- 
liet-near-infrared spectral range. It was found that the coupling 
error must be accounted in the aerosol remote sensing problems 
based on the analysis of the backscattered ultraviolet light. 


Index Terms—Radiative transfer, remote sensing. 


I. INTRODUCTION 


EROSOL remote sensing from space over land is often 

performed assuming that the satellite signal can be decou- 
pled in two components (namely, molecular and aerosol scat- 
tering contributions). Then the lookup tables approach is used to 
retrieve aerosol optical thickness Ta. Such decoupling would be 
an exact procedure as the atmospheric optical thickness Tatm << 
1. Then scattering by molecules and aerosol particles becomes 
independent. However, terrestrial atmosphere is a multiply scat- 
tering medium as far as ultraviolet (UV)—visible—-near-infrared 
(IR) radiative transport is concerned. It means that the total at- 
mospheric reflection function R can be presented as 


R= Rm + Ra + AR (1) 


where Rm is the Rayleigh reflection function defined as the 
atmospheric reflection function for the hypothetic atmosphere 
having only gases and no scattering particles, Ra is the aerosol 
reflection function defined as the atmospheric reflection func- 
tion for the hypothetic atmosphere having only particles and no 
gases. For clarity of considerations, we neglect the ground re- 
flection contribution. However, this could be easily taken into 
account as an additional term in (1). 

The value AR accounts for a coupling between molecular 
and aerosol scattering. Equation (1) allows to present the relative 
error of decoupling process as 


s- AR 
==. 


Calculations show that the value of € is rather small for typical 
atmospheric and satellite observation conditions (usually below 
3% to 5%). However, for aerosol remote sensing applications, 
not coupling error € itself but rather the relative error 6 of the 
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aerosol reflection function determination (in the assumption that 
coupling can be neglected) is of importance. It follows that 
Ra =R} 
§= 4 3 
E, (3) 
where Ra = R— Rm — AR and Rž = R — Rm. We have from 
(3) 


E 


= (4) 


l-e-p 
where p = Rm/R. Clearly, the value of p decreases with A 
and can be neglected for large wavelengths at most sun/satellite 
geometries (e.g., in IR). Then it follows approximately: 6 ~ —e. 
It means that 6 is small if £ is small under assumption p — 0. 
Therefore, 6 and € are of about the same magnitude but of a 
different sign for small p. However, the value of p increases 
toward the UV region and can be in principle considerably larger 
than the coupling error e. In particular, it follows from (4) that 
6 > -easp— 1. 

Let us imagine the situation, when the atmospheric reflection 
function is almost entirely determined by molecular scattering 
processes. Then we have p = 1 — g, where a is a small number. 
Therefore, it follows from (4) 


(5) 


where v = a/ £. It means that 6 — oo as v — 1. Therefore, the 
error of the aerosol reflection function determination neglecting 
coupling could be quite large if R ~ Rm. The coupling errors 
themselves could be very small in this case. Then we need to 
detect a small signal on the much higher level background con- 
tribution. This explains the increase in the error for this case. 

The task of this paper is to study errors £ and 6 for realistic 
atmospheric and illumination/observation conditions. This will 
allow us to understand if one can neglect the coupling contri- 
bution for aerosol remote sensing purposes (e.g., for the deter- 
mination of the aerosol optical thickness from space). The ap- 
proximation assuming AR = 0 looks quite attractive taking 
into account that the gaseous composition of atmosphere is very 
stable. This means that Rm can be calculated quite accurately 
in advance for a given sun/satellite geometry using available in- 
formation on the spectral molecular optical thickness [2]. Then 
the value Ra is approximated as R — Rm, where FR is the mea- 
sured reflection function. Lookup tables of precalculated aerosol 
reflection functions can be searched then to find the correspon- 
dent aerosol optical thickness [12]. 

Similar studies have been performed by Antonie and Morel 
[1] using Monte Carlo calculations at wavelengths larger than 
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TABLE I 
AEROSOL VERTICAL MODEL USED IN CALCULATIONS. OPTICAL MODELS 
FOR EACH AEROSOL TYPE WERE TAKEN FROM HESS et al. [3] AND 
SCIATRAN DATABASE [7]. THE SATELLITE SIGNAL IS MOSTLY 
INFLUENCED BY THE LOWER AEROSOL LAYER 


The position of the aerosol layer Type of aerosol 


0-2km urban 
2-10km average continental 
10-30km sulfate 
30-100km meteoric 
T 
10.0 F — Aerosol 500 nm 4 
a O Aerosol 350 nm 4 
sooo: Rayleigh 500 nm 4 
a 
2 
© 
E 
a 
p oa . 
Do 1.0F 
Ea H 
a H 
0.1 , E E E A apn RS Ai =f 
0) 50 100 150 
scattering angle, degrees 
Fig. 1. Truncated aerosol phase function (50 Legendre moments are used) at 


the wavelengths 350 and 500 nm in the lower layer. The dashed-—dotted line 
gives the Rayleigh phase function. 


445 nm. They concluded that the coupling between aerosols and 
molecules contributes significantly to the path reflectance. We 
extend their study to the case of UV light backscattering. 


II. ATMOSPHERIC AND RADIATIVE TRANSFER MODELS 


In our radiative transfer simulations of reflection functions 
R, Ra, Rm (and, therefore, also 6,¢,1”) we use the following 
atmospheric model. The aerosol optical properties are defined 
in four layers as shown in Table I. The aerosol singl-scattering 
albedos, extinction coefficients, and phase functions have 
been taken from the optical properties of aerosols and clouds 
(OPAC) atmospheric model [3]. The truncated phase functions 
for aerosol scattering in the lower atmospheric layer are given 
in Fig. 1. The truncation is performed according to standard 
techniques [10]. The use of truncated phase functions speeds 
up radiative transfer calculations considerably. 

Aerosol phase functions differ from the correspondent 
Rayleigh phase function also given in Fig. 1. In particular, 
we see that the Rayleigh scattering is much more effective 
in the backscattering region as compared to the aerosol scat- 
tering. This means that purely aerosol and molecular scattering 
media having the same optical thickness will produce different 
backscattering signals. Molecular scattering will be by order 
of magnitude stronger in the single-scattering approximation. 
The spectral single-scattering albedo for the most contributing 
to the reflected light layer 1 (see Table I) is shown in Fig. 2. 
It differs considerably from 1 because of assumed urban 
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Fig. 2. Spectral single-scattering albedo for the lowest atmospheric layer. 


TABLE II 
RAYLEIGH OPTICAL THICKNESS AT SEVERAL WAVELENGTHS 
USED IN CALCULATIONS [2] 


A,nm | 340 370 412 440 550 650 758 865 


T, 0.714 |0.500 |0.319 |0.243 |0.097 |0.049 | 0.026 | 0.016 


soot-polluted aerosol model used in calculations. Because the 
OPAC model is exclusively based on Mie theory, possible 
effects of nonspherisity of aerosol particles are neglected. The 
aerosol optical thickness 7a (550 nm) was varied in the range 
0.05-1.0. The Rayleigh optical thickness have been calculated 
assuming the 1962 US standard atmosphere model with values 
of Rayleigh optical thickness 7, given in Table II. Gaseous 
absorption was neglected. Calculations were performed at 
wavelengths 340, 370, 412, 440, 550, 650, 758, and 865 nm 
using the SCIATRAN 2.0 radiative transfer model developed at 
Bremen University [7]. SCIATRAN 2.0 is based on the discrete 
ordinate technique to solve the scalar radiative transfer equation 
for a plane-parallel turbid layer with account for the vertically 
varying molecular, aerosol, and cloud light scattering. The 
program also incorporates databases of atmospheric air local 
optical characteristics, including gaseous absorption cross sec- 
tions and concentrations profiles. Further details are given by 
Rozanov et al. [7]. The code and correspondent guide are freely 
available at http://www.iup.physik.uni-bremen.de/sciatran. 


IHI. RESULTS OF NUMERICAL CALCULATIONS 


We show the aerosol reflection function R, (solid line) cal- 
culated with account for coupling as compared to the aerosol re- 
flection function R* (dashed line), calculated in the assumption 
AR = 0(no coupling) in Fig. 3. Calculations were performed at 
the spectral points as indicated in Table II at the aerosol optical 
thickness 7,(550 nm) = 0.3, the satellite zenith angle 9 = 30°, 
and the solar zenith angle 99 = 60° for azimuths ¢ equal to 
0°, 90°, and 180°. It follows that R* (no coupling) differs from 
Ra considerably in visible and UV. In particular, the spectral 
aerosol reflection function derived in the assumption that there 
is no coupling takes a maximum around 440 nm and then starts 
to decrease. Such a behavior can be interpreted as an increased 
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Fig. 3. Retrieved aerosol reflection function. 


100 (= 9 aie a aa 
L 8 = 30° ¢ = 90° 


aerosol reflection function error, % 


solar zenith angle, degrees 


Fig.4. Dependence of the aerosol reflection function error 6 on the solar zenith 
angle ĝo. 


aerosol absorption in UV. This is in contradiction with the model 
used (see Fig. 2) and also with experimental results (e.g., see 
[9] and [12]). So clearly the coupling between molecular and 
aerosol scattering cannot be neglected in the satellite retrievals 
algorithms using spectral range below 440 nm (and even below 
500 nm depending on the geometry; see Fig. 3). This is an un- 
fortunate situation taking into account that the reflective power 
of most natural surfaces is much smaller in UV as compared to 
visible. So the UV spectral range is very attractive for aerosol 
remote sensing as compared to retrievals in visible due to the 
minimization of the generally unknown surface influence error 
[5], [11], [13]. 

The angular dependence of the aerosol reflection function 
error 6 is shown in Fig. 4 at @ = 30°, d = 90°, 7,(550 nm) = 
0.3 and wavelengths 340, 370, 412, 440, 550, and 650 nm. 
It follows that the error increases considerably with the solar 
angle and the inverse wavelength. It is around 30% at the wave- 
length 340 nm for the solar zenith angles below 40°. However, 
it sharply increases for low sun producing the positive error 
of about 90%. It means that the retrieved value of the aerosol 
reflectance R* under the assumption of no coupling will be 
too low. This will lead to a considerable underestimation of 
the aerosol optical thickness. The error is smaller at the wave- 
length 412 nm and solar angles smaller than 60°. It is below 
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Fig. 5. Same as in Fig. 4 except for the coupling error £. 


10% and sun-position dependent then. The error is negative and 
small ( typically, below 5%) at wavelengths larger 550 nm and 
solar zenith angles usually used in the aerosol remote sensing 
(o < 70°). Note that aerosol retrievals over oceans are often 
performed in IR channels (e.g., 865 nm) due to the small con- 
tribution of light reflected from ocean in this case (outside of 
the glint region). It is fortunate that the coupling error may 
be neglected for maritime aerosol remote sensing performed at 
A > 700 nm. However, even in this most favorable case, the 
error can reach 15% at the solar angle 80° and the wavelength 
650 nm. This could affect the quality of retrievals in polar re- 
gions, where sun is low. The errors are negligible for wave- 
lengths 758 nm, 865 nm due to low molecular scattering con- 
tribution in this case (see Table II). 

Note that the coupling error ¢ itself is not so large in the cases 
considered (see Fig. 5). The increase of errors toward lower sun 
positions is best explained by the fact that light needs to travel 
larger distances in the atmosphere in this case. This rises up the 
molecular scattering contribution. Indeed, the molecular optical 
depth differs considerably along vertical and along slant paths 
in the atmosphere. This is the reason why e€ is below 1% at the 
solar zenith angle 0° but around —8% at 69 = 80° (see Fig. 5). 
The coupling error is mostly negative suggesting that AR < 0 
in most of cases and, therefore, if one neglects coupling errors, 
too small aerosol reflectances are obtained and, therefore, the re- 
trieved aerosol optical thickness is underestimated. This is con- 
sistent with Fig. 4, where 6 > 0 for most of cases. Also it follows 
from Figs. 4 and 5: 6 ~ —e at small € as it was suggested above. 

Results considered above were calculated for the observation 
angle 30° and the aerosol optical thickness 7,(550 nm) = 0.3. 
We have decided to check how errors depend on the obser- 
vation geometry and aerosol optical thickness. Therefore, we 
have repeated our calculations for other observation geometries, 
and values of Ta. Selected results of this study are shown in 
Figs. 6-8. Results given in Fig. 6 have been calculated at the 
same conditions as in Fig. 4, but the observation angle was 
changed. We see that our conclusion in respect to importance 
of coupling term AR is not influenced very much by the obser- 
vation geometry change. So they are of general nature. 

The influence of the aerosol load is illustrated in Figs. 7 and 
8, where results for 7,(550 nm) = 0.1 and 7,(550 nm) = 1.0 
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Fig. 6. Same as in Fig. 4 except at @ = 0°. 
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Fig. 7. Same as in Fig. 6 except at 7,(550 nm) = 0.1. 
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Fig. 8. Same as in Fig. 6 except at 7,(550 nm) = 1.0. 


are given. Calculations were done at the same conditions 
as in Fig. 6. Again we see that the coupling term cannot be 
neglected in UV. The value of 6 varies in the range 20% to 
90% at A = 340 nm depending on 7, and 69. Dependencies 
on ĝo are monotonous at 7,(550nm) = 1.0 because of greater 
influence of multiple scattering in this case. Errors are generally 
larger then. They are around 20% at the wavelength 412 nm and 


49 < 60°. This will produce the underestimation of the retrieved 
aerosol optical thickness 7,(412 nm) if the coupling is not 
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properly accounted for. The underestimation is even larger at 
0o > 60°, which makes coupling terms even more important 
in this case. This conclusion is of special importance for arctic 
aerosol studies although situations with 7,(550 nm) = 1.0 are 
not frequent there. Fig. 7 shows that the coupling term is still 
of importance even at 7,(550 nm) = 0.1 for the solar angles 
larger than 70° and A < 440 nm. Oscillations on Fig. 7 are 
due to the influence of correspondent oscillations in the phase 
function (see Fig. 1). Note that the scattering angle at the nadir 
observation is equal to 7 — Oo. 

Fig. 9 is similar to Fig. 2 but for another observation con- 
ditions and several aerosol optical thicknesses. It confirms that 
the assumption AR = 0 leads to unphysical maximum in the 
retrieved aerosol reflectance positioned at A ~ 440 nm and gen- 
erally smaller values of the aerosol reflection function as com- 
pared to true (solid lines in Fig. 9) ones. 

Our results are obtained using the scalar approximation of the 
radiative transfer theory. We also estimated the errors related 
to the vector character of light not accounted for in our model. 
For this we have used the vector code SCIAPOL developed by 
us [8]. Polarization effects can be neglected at A > 650 nm 
for Rayleigh scattering and for all wavelengths considered for 
aerosol scattering. However, they are of importance for smaller 
wavelengths for Rayleigh scattering and highly solar angle-de- 
pendent [6]. In particular, we found (not shown here) that errors 
in the Rayleigh scattering component due the scalar approxima- 
tion are in the range 0% to 10% depending on the solar angle in 
the interval of solar angles Aĝo € [0°, 50°], A = 340 nm and 
the nadir observation [8]. The error is smaller for larger solar 
angles in the interval Aĝo. 

It means that instead of (1) we actually have 


R = Rm + Ra + AR + AR, (6) 


where AR, is the contribution to the total reflectance due to 
the vector properties of light. It also means that £, 6 should be 
substituted by 


ee AR+ AR, = e* 
= = 

where the star means that the correspondent error is due to fact 
that both coupling and vector nature of light are not accounted 


(7) 
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Fig. 10. Dependence of the scalar approximation error on the optical thickness 
for Rayleigh scattering at no absorption and various solar illumination angles 
and the nadir observation. 


for. The value of 6 is positive (see Fig. 8) at Ta = 1.0. This 
means that both e = —(1—,)6/(1—6) and AR are negative. We 
found that the same is true for AR, at ĝo > 50° (see Fig. 10). 
This will lead to larger absolute values of the pair (6*,e*) as 
compared to the absolute values of € and 6 and at large solar 
angles. 

Note that it follows: AR, > 0 at 09 < 45° (see Fig. 10). 
In this case the absolute values of the pair (e*, 6*) are smaller 
as compared to the absolute values of £ and 6 due to the partial 
cancellation of errors. This will induce important differences if 
different algorithms are used for the same radiance field inter- 
pretation (e.g., between results obtained using aerosol retrieval 
algorithms in scalar and vector modes at AR = 0). The main 
conclusion is that not only the coupling error but also the polar- 
ization effects should be considered in aerosol remote sensing 
applications in UV (e.g., see [5]). 

UV remote sensing of atmosphere over land is favorable as 
compared to visible-NIR techniques. Then generally unknown 
contribution from the surface is small. However, UV aerosol re- 
mote sensing algorithms suffice from the necessity to account 
for the vector nature of light and increased Rayleigh contribu- 
tion. The latter depends on the ground elevation. Therefore, cor- 
respondent UV retrieval techniques become more complex. This 
is the balanced, however, by much higher accuracy of the UV 
aerosol optical thickness product as compared to Tą obtained 
using visible/NIR retrievals techniques, where light reflected 
from ground makes large and generally unknown contribution 
to the satellite signal. 


TV. CONCLUSION 


We have shown that the coupling of molecular-aerosol scat- 
tering is an important issue of the aerosol remote sensing over 
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land, where wavelengths smaller than 440 nm are used to re- 
duce the uncertainty related to the ground reflectance charac- 
teristics. In particular, the retrieved aerosol optical thickness is 
underestimated in this case. Also a wrong absorption-like fea- 
ture is introduced in the retrieved aerosol reflectance and the 
aerosol optical thickness. This may lead to errors in the retrieved 
aerosol single-scattering albedo, which is of importance for cli- 
mate problem. Results were obtained in a number of simplified 
assumptions like the spherical shape of particles, no ground re- 
flection, and the scalar approximation of the radiative transfer. 
However, main results do not change very much even we ac- 
count for these factors. In particular, nonspherical particles in- 
crease lateral scattering. This will lead to the increase of the 
aerosol phase function. This will make this function closer to 
that of Rayleigh scattering thereby strengthening coupling ef- 
fects. We also estimated the errors related to the vector char- 
acter of light not accounted for in our model. These effects can 
be neglected at A > 650 nm. However, they are of importance 
for smaller wavelengths as discussed in the text. 

One of the important conclusions of this work is that for re- 
trieval of aerosol properties from space for wavelengths below 
440 nm, the iterative vector algorithms as those discussed by 
Hosekamp and Landgraf [4] should be preferably used. The use 
of a priori assumptions on the scalar approximation/strength of 
coupling are not needed then. 

We also show the capabilities of SCIATRAN and SCI- 
APOL radiative transfer codes in the solution of this important 
atmospheric optics problem. The SCIAPOL-_1.0 radiative 
transfer code can be downloaded for noncommercial use at 
www.iup.physik.uni-bremen.de/~alexk. The SCIATRAN_2.0 
code is located at www.iup.physik.uni-bremen.de/sciatran. 
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